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Abstract: 

We introduce the computer code Recola for the recursive generation of tree-level and one- 
loop amplitudes in the Standard Model. Tree-level amplitudes are constructed using off-shell 
currents instead of Feynman diagrams as basic building blocks. One-loop amplitudes are rep- 
resented as linear combinations of tensor integrals whose coefficients are calculated similarly to 
the tree-level amplitudes by recursive construction of loop off-shell currents. We introduce a 
novel algorithm for the treatment of colour, assigning a colour structure to each off-shell cur- 
rent which enables us to recursively construct the colour structure of the amplitude efficiently. 
Recola is interfaced with a tensor-integral library and provides complete one-loop Standard 
Model amplitudes including rational terms and counterterms. As a first application we con- 
sider Z + 2 jets production at the LHC and calculate with Recola the next-to-leading-order 
electroweak corrections to the dominant partonic channels. 
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1 Introduction 



The study of the mechanism of electroweak (EW) symmetry breaking and the search for physics 
beyond the Standard Model (SM) is the primary goal of the Large Hadron Collider (LHC) 
and the corresponding experiments. With the discovery of a bosonic resonance with a mass 
of around 125 GeV important progress has been achieved. Still it remains an open question if 
this resonance is the SM Higgs boson and if there are phenomena of new physics at the TeV 
scale. Evidence for a discovery of new particles and the precise determination of their masses and 
couplings on the one hand as well as the establishment of exclusion limits on the other hand are 
achieved by sophisticated experimental analyses, capable of highlighting a small signal on a huge 
background. For the interpretation of the data a precise knowledge of the background is essential. 
This often relies on theoretical descriptions, sometimes also on data-driven estimations where the 
extrapolation to the signal region is based on theoretical distributions. A sound comparison of 
experimental signals with theoretical predictions allowing precise tests of the SM (or of theories 
beyond) requires high precision from both experiment and theory. 

Theoretical predictions at leading order (LO) in perturbation theory are usually insufficient 
to match the experimental precision. At a hadron collider QCD corrections are indispensable, but 
also EW corrections can have an important impact. For instance, for Higgs-boson production in 
vector-boson fusion, EW and QCD corrections are of the same order of magnitude PQ. Moreover, 
the high energies attained by the LHC allow to collect data in phase-space regions where the 
effects of logarithms of EW origin become sizeable. The high centre-of-mass energies available 
at the LHC generate a lot of events with many particles in the final state. Therefore a proper 
theoretical description of LHC physics requires next-to-leading-order (NLO) computations of 
multiparticle processes (with five, six, or more external legs) in the full SM (including EW 
corrections). 

In the past years many groups have concentrated their efforts to make such calculations 
feasible. New techniques have been proposed, mainly for the computation of one-loop virtual 
corrections which are considered the bottleneck of NLO calculations. In the standard approach 
based on Feynman diagrams the major problem is caused by huge algebraic expressions ap- 
pearing in the computation of the virtual amplitudes. The development of techniques based on 
Generalised Unitarity [2H9] allowed a change of perspective. The formal starting point of these 
methods is the general decomposition of one-loop amplitudes as linear combinations of scalar 
integrals, as obtained from the standard Passarino-Veltman reduction [10J. The computation of 
the coefficients of the scalar functions is then reduced to the calculation of tree-level amplitudes 
by means of cutting equations. The simplicity of these relations allowed the automation of NLO 
QCD computations leading to the development of computer programs [1 1H18] which enabled the 
calculations of many QCD processes [T9W36] of the Les Houches priority list [37] . 

More recently, based on ideas of Soper [38] purely numerical methods for the calculation of 
one-loop QCD amplitudes have been put forward [39]. They are based on subtracting the soft, 
collinear, and ultraviolet divergences of one-loop amplitudes and performing the loop integration 
of the remaining finite integrals numerically after suitably deforming the integration contours in 
the complex space. These methods do not rely on Feynman graphs and have been proven to 
work for multi-parton amplitudes [40j . 

The success of the new methods did, however, not supersede the traditional diagrammatic 
approach. The Generalised Unitarity methods still suffer from the numerical instabilities char- 
acteristic of the Passarino-Veltman reduction. These are overcome by computing the amplitude 
in quadruple or multiple precision at critical phase-space points, reducing however the CPU ef- 
ficiency [HUH]. While rescue solutions in this context have been proposed in Refs. (HUE], it is 
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well-known that numerical instabilities can be avoided by constructing the amplitude as a linear 
combination of tensor integrals, which is the natural representation in the Feynman-diagrammatic 
approach. Various groups j!8 |,l45fH8] have developed efficient computational techniques for the 
calculation of the tensor integrals that avoid numerical instabilities. Using these improved meth- 
ods in the Feynman-diagrammatic approach yields more than competitive numerical codes. In 
particular the methods of Refs. |46j.l48j have been successfully applied to the calculation of both 
EW corrections [Tll4"9H52] and QCD corrections [53ti56] , 

Recently the diagrammatic approach has been further boosted by the advent of the Open- 
Loops algorithm [57]. Organising the diagrams into cut-opened topologies, the coefficients of 
the tensor integrals are recursively built with tree- level- like techniques; the efficiency is increased 
by pinch identities that relate higher-point loop diagrams to pre-computed lower-point diagrams. 
Since the method is based on individual topologies, colour factorises and is treated algebraically. 
The present implementation of OpenLoops can handle NLO QCD corrections to any Standard 
Model process. 

An interesting hybrid method, proposed by Andreas van Hameren in Ref. [58] for evaluating 
one- loop gluonic amplitudes, combines stable and universal tensor-reduction methods for loop 
integrals with a basic result of Generalised Unitarity, namely the reduction of the computation 
of a one- loop amplitude to that of a set of tree-level amplitudes. The technique relies on the 
representation of the amplitude in terms of tensor integrals, whose coefficients are computed 
recursively [59] without resorting to Feynman diagrams at any stage. 

In this paper we present Recola, a generator of SM one-loop (and tree) amplitudes. It is 
based on an algorithm which implements recursion relations for the computation of the coeffi- 
cients of the tensor integrals. The goal is to combine the efficiency of the numerically stable 
tensor-integral reduction with the automation made possible by a completely recursive non- 
diagrammatic approach. 

As a first application of Recola we have calculated EW corrections to pp — > Z + 2 jets. 
Due to its large cross section and similar signatures this process provides a major background for 
Higgs-boson production in vector-boson fusion kinematics. The dominant NLO QCD corrections 
of @{a^a) have been investigated in Refs. [601EI], while a subset of EW Z + 2 jets production 
has been studied at NLO QCD in Ref. [62]. Here, we calculate the EW corrections of 0(a^a 2 ) 
to the dominant partonic channels for Z + 2 jets production. 

This paper is organised as follows: in Section [2] we present the algorithm for the construc- 
tion of amplitudes, introducing first our implementation of recursion relations at tree level 
(Section 12. ip and discussing in Section [2.21 its generalisation to one-loop amplitudes. After some 
remarks on the computation of the rational terms and counterterms (Section I2.3[) . we describe in 
Section [2.41 a new algorithm for the treatment of colour. The setup of our calculation for Z + 2 jets 
production is detailed in Section I3.1[ and numerical results are discussed in Section 13.21 Finally, 
Section [4] contains our conclusions. 

2 RECOLA: REcursive Computation of One-Loop Amplitudes 

Recola is a code written in Fortran90 for the computation of tree and one-loop scattering 
amplitudes in the SM, based on recursion relations. 

2.1 Tree-level recursion relations 

The tree-level recursive algorithm is inspired by the Dyson-Schwinger equations J63H65J and 
follows closely the strategy of Helac [56H5S], using off-shell currents as basic building blocks. 
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Let us consider a process with E external legs, and select a particle P of the modeQ together 
with a sub-set of n external legs. The off-shell current w(P,{n}) is then defined as the sum of 
all Feynman sub-graphs which generate P combining the selected n external particles: 

w(P,{n}) = n \J^y-* ■ 

Here the shaded bubble pictorially represents all possible sub-graphs, and the dot indicates the 
off-shell part of the current. If according to the Feynman rules of the theory a sub-set of n external 
particles cannot generate P, the corresponding current vanishes. For n > 1 the propagator of P 
is included in the definition and the current is called "internal" . A current generated by only one 
(n = 1) external particle P 1 is called "external"; if P = P', it coincides with the wave function 
of the particle, otherwise it vanishes. 

In a theory with tri- and quadri-linear couplings only, the internal off-shell currents can be 
constructed using recursively the Dyson-Schwinger equations: 




i+j=n 

E E :~>^ E E .'OO^:- • ( 2 - 2 ) 



Each term of the sums represents a "branch" , which is obtained by multiplication of the gen- 
erating currents w(Pi,{i}), w(Pj,{j}) [and w(Pk, {k}) for the second term of (12. 2p ] with the 
interaction vertex, marked by the small box, and the propagator of P, marked by the thick 
line. Branches have to be built for all possible generating currents formed by sub-sets {i}, {j} 
(and {k}) of the n external particles such that % + j (+k) = n. The recursive construction of 
the amplitude can be efficiently implemented using lower-multiplicity off-shell currents as seeds 
for the numerical evaluation of higher-multiplicity ones. First, all possible currents with two 
external legs (n = 2) are calculated combining through a tri-linear coupling pairs of external 
currents only. Next the currents with three external particles (n = 3) are generated summing 
the branches with three external currents linked through a quadri-linear coupling and those with 
one external current and one of the calculated internal currents with two external legs, linked 
through a tri-linear coupling. Analogously the currents with n = 4 are then computed combining 
two or three currents with n = 1, n = 2 and n = 3, and so on. As a consequence of the summation 
in (|2.2p . the current w(P,{n}) depends on the particle P and on the set of generating external 
particles {n} but not on the particular way these particles have been combined in order to get 
w(P, {n}). In this respect, working with off-shell currents rather than Feynman diagrams allows 
to avoid recomputing identical sub-graphs contributing to different diagrams, since each current 
is computed just once; furthermore, the summation in ()2.2[) reduces the number of generated 
objects which are passed to the next step of the recursion. 

In order to calculate the amplitude for a tree-level process A — > B we first use crossing 
symmetry to switch to incoming particles and consider the corresponding process A + B — > 0, 
where B is charge conjugate to B. Next we choose one of the external legs, for definiteness 
the .E th leg, to close the construction of the amplitude. Then all currents resulting from the 
other E — 1 external legs are recursively constructed. In the last step of the procedure, i.e. for 



1 Here all particles of the model have to be taken into account, also unphysical ones like would-be Goldstone 
bosons. 



3 



n = E — 1, we require the generated particle P to be the selected E external particle. As a 
consequence this last current is unique, and the amplitude M is obtained multiplying with the 
inverse of its propagator and with the wave function of the selected E th particle (which coincides 
with the E th external current): 

M = E-lty^y. x (propagator)^ 1 x . E . (2.3) 

The recursive evaluation of the currents begins with the external currents (re = 1), which, for 
colourless particles of a given polarisation A and momentum p, are given by the corresponding 
wave functions: 

A-^-* = u x (p), = Am = e x (p), = 1. (2.4) 

p p p p 

The explicit expressions for the spinors u\(p), u\(p), v\(p), v\(p) of fermions and for the po- 
larisation vectors e\(p), e* x (p) of vector bosons have been coded using Ref. [69]. For coloured 
external particles also the information on colour has to be kept, as we explain in Section [2.41 

Once the external currents have been numerically evaluated, the internal ones are built using 
the Feynman rules of the theory. For example, given a pair of external e~ and e + with momenta 
pi and pi and currents Ux 1 (pi) and V\ 2 (p2), one can generate the internal current of a photon 
contracting the two external currents with the QED vertex — ie7^ and the photon propagator 
— ig^u/ipi +P2) 2 - Once the particle P of the off-shell propagator is fixed, in general, several 
branches contribute to the same internal current since different combinations of particles with 
appropriate interaction vertices can generate the same particle P. In these cases the contributions 
can be simply summed up, according to (|2.2p . 

The recursive algorithm is implemented in the code Recola through two steps. 

In the first part, the initialisation phase, the currents are identified by integer numbers which 
contain all the relevant non-dynamical information: the particle content (e _ , e + , 7, etc.), the 
colour information (see Section [2.4p and an integer tag number. The tag number is assigned 
according to a binary notation |66|I70|. In practice, the external currents get a label 2 l_1 , where 
i is the i th external leg, i.e. 1 — y 1, 2—7-2, 3—7-4, . . . , E — > 2 E ~ 1 . The tag number of an internal 
current is obtained by a summation of the integer tags of the external currents contributing to it. 
The binary notation ensures that, given the tag number of any internal current, the contributing 
external legs can be uniquely identified. Moreover, it reflects the basic property that a current 
depends on the generating external particles {re}, but not on the particular way these particles 
have been combined in order to obtain the current. 

The initialisation part of the code builds a skeleton of the amplitude: all needed off-shell 
currents are enumerated and, for each branch, all generating off-shell currents and the generated 
one are identified. This part is run once for all, before giving explicit values to the momenta 
of the external particles, i.e. before performing the phase-space integration in a Monte Carlo 
program. 

The second part of the code, the dynamical production phase, uses the results of the first 
part to actually compute the amplitude for each point of the phase-space. First, the external 
currents are numerically computed; then, branch after branch, all internal currents are recursively 
evaluated according to the skeleton generated in the initialisation phase. Here the code has to 
be as efficient as possible in terms of CPU time because it must run on a large grid of points in 
phase-space; the computation of all non-dynamical quantities in the initialisation phase allows 
to avoid a repetition of those operations which can be performed once independently of the 
particular values of the external momenta. 
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2.2 One-loop recursion relations 



Let us now move to the one-loop case. After summing all contributing Feynman graphs Qi every 
one-loop amplitude can be written as a linear combination of tensor integrals: 



Here the tensor coefficients cL'..^. do not depend on the loop momentum q, which is present 
only in the tensor integrals 



The index j classifies the different tensor integrals needed for the process, the integer Nj equals 
the number of loop propagators, and r, ( < Nj in the 't Hooft -Feynman gauge) is the rank of 
the tensor integral. 

Leaving aside the computation of the tensor integrals, to be performed with the preferred 

(j r ■ N-) 

technique, we focus here on the tensor coefficients c^.. 3 .^ , which for multi-leg processes, due 
to the complexity of the SM, result in long algebraic expressions in the standard Feynman- 
diagram approach. An interesting idea has been proposed in Ref. [58J, where recursive relations 
for the tensor coefficients of gluon amplitudes have been derived for colour-ordered amplitudes 
of purely gluonic processes. We have further developed this approach to deal with the full 
SM. There is a clear topological correspondence between a one-loop diagram with E external 
legs and a tree diagram with E + 2 external legs, obtained after cutting one of the loop lines. 

(j r ■ N) 

After uniquely fixing this correspondence, one can compute the tensor coefficients Cp.l---^. with 
recursion relations similar to those used for tree amplitudes. 

Given a process A — > B at one loop, we consider first the set of all tree processes A + B + 
P + P -> for each particle P of the SM@ However, the set {A + B + P + P -> 0, VP G SM} 
contains more diagrams than the original one-loop process A — > B. This is due to the fact that 
we can cut the loop diagram at any of its loop lines and that we can run along the loop clockwise 
and counterclockwise. Therefore, we have to fix some rules to discard the redundant diagrams. 
To explain these rules, we can work without loss of generality in a theory with a single scalar 
particle (j) with a tri-linear coupling (the generalisation to the presence of a quadri-linear coupling 
is straightforward). In such a theory the set {^4 + B + P + P — > 0, VP € SM} reduces to one 
process, i.e. A + B + <j) + (j) — > 0. We use tag numbers for external and internal currents as 
explained in Section 12.11 and assign tag numbers 2 E and 2 E+l to the currents corresponding to 
the two additional external legs of P and P. These currents are called "external loop currents" , 
while the legs of P and P are called "external loop legs". 

Let us first consider the sets of diagrams with three and four external legs where only external 
particles enter the loop. Marking with a cross the two external loop legs of the trees, we get: 



4 8 16 8 16 8 16 8 16 8 16 8 16 



(2.5) 





(2.6) 



2 Here all unphysical particles, in particular also Faddeev-Popov ghosts, must be included. 
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The diagrams on the right-hand side are obtained by cutting in all possible ways one of the loop 
lines of the diagrams on the left-hand side. The tree diagrams have been drawn in such a way, 
that one can easily identify the original sequence of the loop lines (called "loop flow"), starting 
from the external loop leg with tag number 2 E and ending with the external loop leg with tag 
number 2 E+1 . Two simple rules can avoid the redundant tree diagrams and fix properly the 
correspondence with the loop diagrams: 

1') The external current 1 must be attached to the external loop current with tag number 2 E . 
This rule fixes the starting point of the loop flow and thus reduces the redundancy already 
up to a factor of two, the direction of the loop flow. 

2') The external currents 1, 2 and 4 must be attached to the loop flow in ascending order (the 
other external currents can enter the loop flow everywhere, also between 1, 2 and 4). This 
rule uniquely fixes the direction of the loop flow. 

With these rules the way of cutting each loop diagram of (|2.8[) becomes unique: 

4 8 16 4 



2 2 

1 g 1 4 1 4 1 16 32 8 1 16 32 4 1 16 32 4 




2 42 88 22 42 88 2 



The rules have to be generalised to other classes of diagrams where external legs can combine 
in tree sub-graphs before entering the loop flow. To this end we define an identifier number 
for each current, given by the smallest external tag among those forming its tag number. For 
example, a current with tag number 13 = 1 + 4 + 8, which has been created combining the 
external legs 1, 4 and 8, has identifier 1; a current with tag number 6 = 2 + 4 has identifier 2. For 
external currents the identifier coincides with the tag number. In case of a quadri-linear coupling, 
the two currents entering the loop are represented by a common identifier, the minimum of the 
two identifiers for the individual currents. 

Now the generalisation of the two rules is straightforward: 
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1) The current with identifier 1 must be attached to the external loop current with tag number 
2 E . 



2) The currents with the three smallest identifiers must be attached to the loop flow following 
the ascending order of their identifiers. 

The proper treatment of self-energy insertions deserves particular care. For two identical 
particles flowing in the self-energy loop, the selection rules (actually the first rule alone) reduce 
the number of cut diagrams to one. On the other hand in this case the self-energy diagram 
gets a symmetry factor 1/2 because of Wick's theorem. If the two particles in the loop are 
different, we get two cut diagrams, which however give the same contribution. Therefore, the 
cutting procedure for self-energies reproduces the loop diagrams correctly if the corresponding 
tree contributions are multiplied by a factor 1/2 in all cases. 

In addition to rules 1) and 2), in most renormalisation schemes some classes of diagrams have 
to be discarded, namely tadpoles and self-energy insertions on external legs. Therefore, being 
S = 1 + 2 + 4 + . . . + 2 E ~ l = 2 E - 1 the sum of all external tags of the process, we use the following 
additional rules: 

3) A current with tag number equal to S or equal to S — 2 n , n = 0, 1, . . . E — 1, cannot enter 
the loop flow in a branch with a tri-linear vertex. This eliminates tadpole diagrams and 
those self-energy contributions made of tri-linear vertices which are inserted on external 
legs. 

4) If in a branch with a quadri-linear vertex one of the two currents entering the loop flow is 
external, the sum of their tag numbers cannot be equal to S. This eliminates the self-energy 
contributions involving one quadri-linear vertex inserted on external legs. 

Applying these four rules to the generation of the skeleton for the off-shell currents of the tree 
processes of the set {A — > B + P + P, VP £ SM}, we obtain the proper skeleton for the one- loop 
process A — > B. 

Having reduced the formal generation of the one-loop amplitude to the generation of a set of 
tree-level processes, we can build the "loop off-shell currents" in a similar way as the tree- level 
currents in Section [2~T| in order to obtain the tensor coefficients c^^.l'^. of (|2.5p . At one- loop 
level, the particle which closes the recursive construction is the external loop leg with tag number 
2 E+l . Since all loop lines are virtual lines and retain their propagator, the last step of (|2.3p . 
where the last generated current is multiplied with the wave function of the particle closing the 

recursion, • x 2 E+1 , is performed without multiplication by the inverse propagator resulting 

in 

2 E 



5M = 




(2.10) 



1 2 E ~ l 1 2 s " 1 



The external currents for the first E external legs are defined as in Section [2.1t the external 
currents of the two external loop legs are defined such that the contraction originally contained 
in the loop can be easily reproduced. To this end, we introduce a suitable set of spinors tpi = Ui,V{ 
and polarisation vectors for the cut fermions and vector bosons, 

4 

(tpi)a = $i)a = $ia, with ^(tfi) a(A) = <W> 

i=l 
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6f = 5f, With £ e f< = ^ (2 . n) 



2=1 



where i denotes the "polarisation" , and a, j3 and /i, v are spinor and Lorentz indices, respectively. 
The loops are glued together as: 



1 



scalar s: \ — s -H- 





vector bosons: \ — 5 -H- > t — I x *aaa/sx£ 



fermions: t 3 / c 5 x »-»-xV ; 




«-> > tzS x «^«-x^. (2.12) 



Except for the scalar case, the cutting procedure associates to the one-loop amplitude the sum 
of four tree-level amplitudes with particular spinors/polarisation vectors for the cut particle. 

Having fixed the external currents, we describe how to compute the internal ones. As ex- 
plained in Section 12.14 f° r tree amplitudes these are computed summing up the currents generated 
in branches where the generating currents are multiplied with the Feynman rules for the vertex 
and the propagator of the generated particle. This is valid also at one-loop level for pure tree 
currents built by combining the original external legs 1, . . . , 2 E ~ l . 

The new features of the loop case are connected to the loop off-shell currents involving the 
external loop leg with tag number 2 E carrying a loop momentum q. The external loop current 
with tag number 2 E defines the beginning of the loop flow; all currents with tag number > 2 E 
belong to the loop flow and are called loop currents, while the branches generating them are 
called loop branches. Every internal loop current contains a g-dependence, generated by the 
Feynman rules for the vertex and the propagator. Working in the 't Hooft-Feynman gauge in 
the SM, the g-dependence of (vertex) x (propagator) takes the form 

(vertex) x (propagator) = £5 o, (2-13) 

[q + p) — 

where the linear g-dependence in the numerator results from a fermion propagator, from a cou- 
pling of three vector bosons, or from a coupling between one vector boson and two scalar/ghost 
particles (in all other cases = 0). If the interacting particles are fermions and/or vector 
bosons, the coefficients and b have an additional Dirac and/or Lorentz structure which is not 
made explicit here for simplicity. Denoting by w\{q) the first internal loop current, we have 



(q+Pi) 2 - m\ 



mi) = ,:r:\ 2 s ^ 



where p\ is the sum of the external momenta entering the first loop branch while and 
d\ t o result from a product of the generating currents with the constants a M and b in (|2.13p . 
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If the current w\ corresponds to a fermion/ vector boson, w\ as well as d±\ and dig carry 
an additional spinor/Lorentz index, suppressed in (I2.14|) . Proceeding along the loop flow, the 
second internal loop current is built combining w\(q) with tree currents and with a product 
(vertex) x (propagator) of the form (12. 13H and reads 



4,2^1 ffM2 + + d 2,o 

[(<? + Pi) 2 - m i] [(q + P2) 2 - m\ 



= r,, ^2 ]Un + ^2_^2 V ( 2 - 15 ) 



The recursively constructed I loop current of the loop flow is of the form 

Wl (q) = y\ fX" 1 *^-** (216) 
fc = nLi[(9+^) 2 -m2] 

For the last loop current (with I = Nj = number of loop lines) the momentum p^. is equal to 
the sum of all external momenta and thus vanishes. 

The recursion relations (|2.2p are valid also for the loop currents, but cannot be used to com- 
pute them numerically (unless we give an explicit value to the loop momentum q). However, one 
can define similar relations to compute the set of coefficients {dip, df\, . . . , dfj" w }, using the gen- 
eral form of (|2.13p for the g-dependence of loop branches. In fact in a loop branch with a tri-linear 
vertex, knowing the generating tree current wt and the coefficients {di-\ t o, d 1 ^^ 1 , . . . , d^ 1 j^ 1 } of 
the generating loop current, the coefficients {d/^,^^, • • • ,df^ of the generated loop current 
are given by 

Ho, dg, . . . , dg-«} = vh ({0, 4_ 1>0) dfc hl , dptfLJ aw 

+ {d K „,C 1 , 1 ,..,CiH.0}4 (2.17) 



where we have again omitted the Dirac/Lorentz indices associated with fermionic or vectorial 
currents as in (|2.14p . For quadri-linear vertices, the situation is even simpler because in this case 
a' 4 = and <[27T7|) simplifies to 

{di,o,<i > ... J <r W > = w tl w t Adi-i,o,d^ lv ...,d^ 1 Ji- 1 1 ,0}b, (2.18) 

where and wt 2 are the two generating tree-level currents. The expressions (|2.17p and (12.181) 
have to be used at each loop branch. The generated coefficients dfi ' , k = 0, . . . , I are in general 
not symmetric under the exchange of their indices fi\ ■ ■ ■ but, being implicitly multiplied by 
the symmetric product q^ ■ ■ ■ q^ k , they can be symmetrised at each step. In this way the number 
of independent coefficients df I is decreased, leading to a reduction of operations in subsequent 
steps of the recursion. 

The recursion relations for loop branches allow us to compute the coefficients of the loop cur- 
rents, but we are not allowed to sum them unless their denominators are equal. From (|2.16p one 
can see that the denominators are products of propagators and are determined by a sequence of 
off-set momenta {pi, . . . ,pi} and masses {mi, . . . , mi}. Therefore, while tree currents are defined 
by the tag number, the particle content, and the colour information, the loop currents need an 
additional parameter, called sequence number, containing the information on {pi, ■ ■ ■ ,pi} and 
{mi, . . . , mi}. in this way, loop currents with different g-dependent denominators are distin- 
guished, and contributions from loop branches with the same denominators can be summed as 
for tree branches. 
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The introduction of the sequence number spoils the uniqueness of the last current. Given the 
cut particle and its polarisation [the index i in (|2.1ip ], the coefficients {d^ 0) dfj. 1 , . . . , dQ. } 

3 1 3 ' ' J ' 3 

of the last current with sequence number n s give a contribution to the tensor coefficients of (|2.5p 
with j = n s . The index j, describing the class of the tensor integrals, is then identified with the 
sequence number of the last currents. At this level, contributions from different polarisations of 
the cut particle can be summed up. For different cut particles one gets in general contributions 
to different classes of tensor integrals. However, since not the particles but only their masses 
and momenta enter the sequence number, also contributions to the same tensor integral classes 
appear and are combined. 

Also at one-loop level the code is divided into an initialisation and a production phase. In 
the initialisation phase, the skeleton of the branches is generated and all quantities which do not 
depend on the momenta are fixed. In particular, the sequence numbers of the last currents allow 
already at this step to determine the list of needed tensor integrals. Therefore the computation 
of the tensor integrals can be done independently of the one of the tensor coefficients (in the 
production phase of the code). 

2.3 Rational terms and renormalisation 

In dimensional regularisation the calculation of Feynman amplitudes is performed in D = 4 — 2e 
space-time dimensions, and the result is arranged as a power series in e. In this way, UV 
divergences of the loop integrals manifest themselves as 1 /e poles to be subtracted upon renor- 
malisation. For consistency, all space-time related objects entering the amplitude have to be 
promoted to their D-dimensional generalisations and all manipulations have to be performed in 
D dimensions. Otherwise terms of order 0(e) are missed which, combined with the 1/e pole of 
the loop integral, give finite contributions to the amplitude. Because these are rational functions 
of the kinematical invariants, they are conventionally called rational terms. A rational term is 
dubbed i?i-term if it results from the e-dependence of the denominators of the loop integrals 
and it is called ii^-term if it is generated by a 0(e) term in the numerator of the Feynman 
amplitude }7T] . 

We assume that the tensor integrals, taken as input by Recola, contain the i?i-terms, either 
by keeping the denominators of the loop propagators D-dimensional in the calculation or by 
explicitly adding these terms. Note, however, that even if the calculation of the tensor integrals 
is performed in D dimensions, they enter Recola as numerical four-dimensional tensors. The 
numerical construction of the tensor coefficients, on the other hand, works strictly in four space- 
time dimensions, so that the i?2-terms are not taken into account automatically. These terms 
can, however, be easily computed using effective Feynman rules which have been implemented 
in our code using the results of Refs. [71H74j l^l. The insertion of the effective Feynman rules for 
vertices and propagators is performed in the tree-level amplitude generator, taking care that only 
one of the vertices results from a rational term. 

Renormalisation is performed via counterterms based on the conventions of Ref. |75j . In 
analogy with the effective Feynman rules for the rational terms, the insertion of counterterms 
takes place in the tree- level amplitude generator. Presently, counterterms are fixed following the 
complex-mass scheme of Refs. [761177] and the results of Ref. [75] for all parameters of the SM. The 
strong coupling constant is renormalised in the MS-scheme at a general scale Q for contributions 
coming from gluons and light quarks, while the top-quark contribution is subtracted at zero 
momentum. 

3 We thank R. Pittau for clarifications concerning Refs. |73II74| . 
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2.4 Treatment of colour 



In the computation of the currents an important aspect is the treatment of colour, which does 
not factorise in the recursive construction (contrary to the diagrammatic approach). One could 
obtain the factorisation of colour by splitting the amplitude in a sum of colour-ordered ampli- 
tudes, as done in Ref. [78j . This, however, would increase the number of amplitudes to compute 
and would become complicated in the full SM. Alternatively, one could compute colour-dressed 
amplitudes, as for instance in Ref. [66], where off-shell currents would carry explicit colour in- 
dices and would have to be computed for each index separately, slowing down the calculation 
considerably. Although the number of colour-dressed amplitudes to be computed can be de- 
creased by a Monte Carlo sampling over colour configurations [67|, the number of operations at 
intermediate steps remains large. In order to optimise the colour treatment further, we devel- 
oped an alternative approach based on "structure-dressed" amplitudes, where each current gets 
an explicit colour structure. 

This is easily achieved working in the colour-flow representation of the 1/N C expansion [79] . 
introduced in Refs. [80j[81] for perturbative QCD computations, where the conventional 8 gluon 
fields are replaced by a 3x3 matrix (A^fj = 4= A^X 11 )^ with the trace condition J2i (A*)i = 0- 
Quarks and antiquarks maintain the usual colour index i = 1,2,3, while gluons get a pair of 
indices i,j = 1,2,3; the Gell-Mann matrices A a and the structure constants in the Feynman 
rules are then substituted by combinations of Kronecker 5s. The propagators read 
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while the vertices become 
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In all Feynman rules the colour part is described by products of Kronecker 5s, and therefore in 
the colour-flow representation the colour structure of the amplitude can be simply obtained as 
a linear combination of all possible products of Kronecker 5s carrying the colour indices of the 
external particles. For a process with k external gluons and m external quark-antiquark pairs 
the amplitude takes the simple form: 



tai,...,a n 



P(l,..,n) 



n 



k + m, 



(2.21) 



where in general all n! permutations P(l, . . . ,n) of the indices /3±, . . . , f3 n have to be considered. 

In a framework based on colour-dressed amplitudes, the colour indices of the external par- 
ticles would be fixed and at each branch, given the colour of the generating currents, all colour 
configurations (3 for quarks or antiquarks, 9 for gluons) for the generated current would be 
computed. Many of them are zero, and the others differ just by simple factors. Since in this 
approach one would define and compute unnecessarily many currents, we decided to follow a 
different strategy. 

Instead of assigning an explicit colour to the currents, we assign them a "colour structure", 
which is a product of Kronecker 5s. In order to understand how these structures look like, let us 
first consider the external currents for quarks, antiquarks, and gluons: 



a- 



a- 



vx(p)Sf, 
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(2.22) 



where a and j3 are the colour indices of the external particles while i and j are "open" colour 
indices, which, during the recursive construction of internal currents, are contracted with the 
indices of the Feynman rules of ([2.21|) . In the recursive procedure these contractions generate 
products of 5s: some of them carry indices of external particles only, as in (|2.2ip . and some 
others involve the open indices of the generated current. For example, the combination of an 
external quark with colour structure S 1 } with an external gluon with colour structure 5 % ^ 5^ 
produces, according to the Feynman rules of (I2.21|) . a quark with two possible colour structures: 
5g 2 5 l and 5g 2 5 l . In both structures the first 5 carries just external indices a±, fa, while 



the second one contains also the open index % of the generated quark current. 

The resulting colour structures for the off-shell currents are in complete correspondence to the 
colour structure in (|2.21j) for the full amplitude. The indices of the ^-structures of a particular 
off-shell current are given by the colour indices a, , j3j of the external particles generating the 
current and by potential open indices for the generated particle: no open indices for a colour- 
neutral particle, one for a quark/antiquark, two for a gluon. Therefore, in general, the colour 
structure of a gluon current, obtained from k external gluons and n — k external quark-antiquark 
pairs, takes the form 




~!5i 5f n 



(2.23) 
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where permutations P(l, ■ ■ ■ ,n + 1) of the indices /3±, . . . ,/3 n ,j on the right-hand side correspond 
to different currents. For the colour structure of a quark (antiquark) current, obtained from k 
external gluons, n — k — 1 external quark-antiquark pairs and an additional quark (antiquark) 
we have 



3 S^Zl 6^,(2.24) 



where again permutations P(l, . . . ,n) of . . . , (3 n in the first case and of /3\, . . . , /3 n -\,j in the 
second correspond to different currents. 

We can easily distinguish two parts in the colour structure: the "open part" , which is always 
present in coloured currents, contains one (for quarks and antiquarks) or two (for gluons) 6s with 
open colour indices i and/or j, while the "saturated part", which is absent in external currents, 
is a product of 5s with only external colour indices. Only the open parts of the structures play 
an active role in the combination of currents, while the saturated parts of the generating currents 
simply multiply to give the saturated part of the generated current. In our code these two parts 
are represented by integer numbers, based on a binary notation. 

Moving from colour-dressed to structure-dressed currents reduces already the number of 
currents. For example, there are 9 colour-dressed currents for a gluon generated from the currents 
of a quark and an antiquark (although many of them vanish), while we have just one structure- 
dressed current. A further optimisation can be achieved introducing a colour label and giving the 
same colour label to currents differing just by a colour factor (due to subsequent multiplication 
of different colour coefficients). In the example considered above of an external quark combined 
with an external gluon, the coefficients of the two possible colour structures <5^<5^ 2 and 5^ ^ 
differ only by a factor —1/N C . For structure-dressed currents this label is easily introduced, 
together with the corresponding colour factors, already in the initialisation phase and can be 
used in the production phase of the code to compute just one of the currents with the same 
colour label. 




3 Electroweak corrections to Z + 2 jets production at the LHC 

As a first example for the application of the code Recola, we consider the EW corrections to 
the dominant partonic channels contributing to the process pp — > Z + 2 jets. 

3.1 Details of the calculation 

At leading order (LO) in perturbation theory, the production of a Z boson at the LHC in 
association with a pair of hard jets is governed by the partonic subprocesses 

Qig-tQigZ, (3-1) 
Qi Qi -» Qj Qj Z, qu qj = u, c, d, s, b, (3.2) 

and their crossing-related counterparts. Since we neglect flavour mixing as well as finite-mass 
effects for the light quarks, the LO amplitudes do not depend on the quark generation, and the 
contributions of the various generations to the cross section differ only by their parton lumi- 
nosities. While the mixed quark-gluon (gluonic) channels (|3.ip contribute to the cross section 
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Figure 1: From left to right: Sample tree diagrams for the QCD contributions to qig — > qigZ 
and to qi q~i — > qj q~j Z, and the EW contributions to qi q~j — > qi q~j Z. 

exclusively at order 0(aa 2 ), the four-quark channels (|3.2p develop LO diagrams of strong as 
well as of EW nature leading to contributions of order 0{aa 2 ), 0(a 2 a s ), and (D(a 3 ) to the cross 
section. Representative Feynman diagrams are shown in Fig. [TJ If standard experimental accep- 
tance cuts are applied (see Section 13.2.11 for the specification of our cuts) , the mixed quark-gluon 
channels clearly dominate over the four-quark channels, with the subprocesses ug — > ugZ and 
dg — > dgZ contributing ~ 70% and the complete class (|3.ip of partonic subprocesses contributing 
~ 80% to the total cross section^ Therefore as a first step towards a complete NLO calculation 
of EW effects in Z + 2 jets we here calculate EW corrections to the gluonic channels (|3.ip . 

3.1.1 General setup 

In our calculation we describe potentially resonant Z-boson propagators appearing in loop dia- 
grams (see Fig. [2] left for a sample diagram) by attributing a complex mass 

ix\ = Ml - iAf z r z (3.3) 

to internal Z bosons. To this end, we consistently use the complex-mass scheme |76 p 77 t l82] where 
fj, 2 ^ and /i| are defined as the poles of the W- and Z-boson propagators in the complex plane. On 
the other hand, the external Z boson is treated as a stable final-state particle with its invariant 
mass being fixed to M z , where M| is the real part of the complex pole of the Z-boson propagator. 
The pole values My and Ty (V = W, Z) for the mass and width of the W and Z boson are related 
to the on-shell results My S and Ty S obtained from the LEP and Tevatron experiments by [83] 

m v = m^/^i + (ros/ M os)2 ( r y = ros /v /i + (ros/Mos) 2 . (34) 

For the definition of the electromagnetic coupling constant a we adopt the scheme, i.e. 
we fix the value of a via its tree- level relation with the Fermi constant G^: 

V^M^f M^\ 

Compared to the Thomson-limit definition of a, the definition of a^^ in the G^-scheme incorpo- 
rates effects of the renormalisation-group running from the scale Q 2 = to the scale Q 2 = M^. 
In addition NLO corrections involving logarithms of light quark masses are avoided as such 
contributions do not enter the muon decay. 

4 In a scenario where vector-boson fusion kinematics is imposed, the dominance of the gluonic channels (13.11 1 
shrinks. Requiring, for instance, two tagging jets in opposite hemispheres with rapidity difference ly^ — yj 2 \ > 4, 
the gluonic channels contribute only ~ 65%. 
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Figure 2: Left: Box diagram involving a potentially resonant Z-boson propagator. Right: Pen- 
tagon diagram involving a 5-point tensor integral of rank 4. 

3.1.2 Virtual corrections 

The virtual corrections involve 0(300) diagrams per partonic channel, including 20 pentagon and 
71 box contributions H The most complicated topologies are given by pentagons involving 5-point 
functions of rank up to r = 4 (see Fig. [2] right for a sample diagram). The virtual amplitude 
is calculated using the 't Hooft-Feynman gauge. The calculation of the tensor integrals is per- 
formed employing recursive numerical reduction to scalar integrals based on Refs. [4*6 l l4^ 1[84Tl86] . 
Numerical instabilities from small Gram determinants are avoided by resorting to various ex- 
pansion algorithms for the problematic momentum configurations |48| . Both, in the case of UV 
divergences as well as in the case of infrared (IR) divergences, dimensional regularisation is ap- 
plied to extract the corresponding singularities. The EW sector of the SM is renormalised using 
an on-shell prescription for the W- and Z-boson masses in the framework of the complex-mass 
scheme |77j . As the coupling is derived from M\y, M% and G^, its counterterm inherits a 
correction term Ar from the weak corrections to muon decay |87j . 

For virtual NLO contributions the finite top-quark mass affects partonic channels involving 
external bottom quarks in a different way than channels with external quarks of the first two 
generations. While the top-quark mass is properly taken into account in closed fermion loops, 
finite top-quark-mass effects constrained to diagrams with external bottom quarks are neglected 
(bg — > ( bgZ and bb — > ggZ are suppressed by the bottom PDFs, gg — > bbZ contributes about 
1% at LO). 

3.1.3 Real corrections 

The EW real corrections to the subprocess (|3. 1[) are induced by photon Bremsstrahlung and 
given by 



Emission of a soft or a collinear photon from an external quark leads to IR divergences which 
are regularised dimensionally. If an IR-save event definition is used, the final-state singularities 
cancel with corresponding IR poles from the virtual corrections. For the initial-state singularities 
this cancellation is incomplete but the remnant can be absorbed into a redefinition of the quark 
distribution function. Technically we make use of the Catani-Seymour dipole formalism as 
formulated in Ref . |88j , which we transferred in a straightforward way to the case of dimensionally 
regularised photon emission. 

In addition to the singularities from soft and collinear photon emission we face a further 
source of IR divergences originating from a soft final-state gluon (see Ref. [50]). Isolated soft 

5 While Recola does not use Feynman diagrams, we give these numbers as a measure of the complexity of the 




(3.6) 



process. 
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gluons do not pose any problem as they do not pass our selection cuts because the requirement of 
two hard jets is not fulfilled. However, in IR-safe observables quarks, and thus all QCD partons, 
have to be recombined with photons if they are sufficiently collinear. Thus, if a soft gluon is 
collinear to a photon, it still passes the selection cuts if recombined with the collinear photon, 
giving rise to a soft-gluon divergence that would be cancelled by the virtual QCD corrections to 
Z + l jet + 7 production. Following Refs. |50U89| we eliminate this singularity by discarding events 
which contain a jet consisting of a hard photon recombined with a soft parton a (a = qi,q~i,g) 
taking the photon-jet energy fraction z 7 = E^/(E y + E a ) as a discriminator. Photonic jets with 
z 7 above a critical value z^ ut are attributed to the process pp — > Z + 1 jet + 7 and therefore they 
are excluded. However, this event definition is still not IR-save because the application of the z 7 - 
cut to recombined quark-photon jets spoils the cancellation of final-state collinear singularities 
with the IR-divergences from the virtual corrections. This is cured by absorbing the left-over 
singularities into the measured quark-photon fragmentation function |90|I91]. 

3.1.4 Implementation 

In order to ensure correctness of our results, and in particular of the calculation with Recola, 
we have performed two independent calculations which we find to be in mutual agreement. While 
the first one applies the technique of recursive amplitude generation as described in Section [21 
the second one relies on the conventional Feynman-diagrammatic approach. 

In the first calculation the amplitude generator Recola provides the Feynman amplitudes. 
For the evaluation of the tensor integrals Recola is interfaced with the Fortran library Col- 
lier [92J. To this end Collier has been extended by an efficient algorithm for building up 
the tensor integrals from the recursively calculated Lorentz-invariant coefficient functions. The 
phase-space integration is performed by means of a generic in-house Monte-Carlo generator [93] 
following the multi-channel sampling approach. 

The second code uses FeynArts 3.2 [MIES] and FormCalc 3.1 [96] for the generation 
and simplification of the Feynman amplitudes. For the numerical evaluation the amplitudes 
are translated into the Weyl-van der Waerden formalism [97] using the program Pole [98] . 
The tensor integrals are again evaluated by Collier which by itself provides two independent 
implementations of all its building blocks. Finally, the phase-space integration is performed with 
the multi-channel generator Lusifer 99 . 

3.2 Numerical results 

3.2.1 Input parameters and selection cuts 

We use the following set of input parameters |100] , 

Gf, = 1.1663787 x lO^GeV" 2 , 

M$ s = 80.385 GeV, r^ s = 2.085 GeV, 

Mg s = 91.1876 GeV, Tg s = 2.4952 GeV, 

M H = 125 GeV, m t = 173.2 GeV (3.7) 

with the value for the top-quark mass taken from Ref. [101j . The pole masses A%,z and widths 
I\y,z entering our calculation are obtained from the stated on-shell values M^ s z according to 
(|3.4j) . The electromagnetic coupling constant ac M is determined from G^ via (|3.5p . The CKM ma- 
trix only appears in loop amplitudes and is set to unity. 

For the prediction of the hadronic pp — > Z + 2 jets cross section the partonic cross sections 
have to be convoluted with the corresponding parton distribution functions (PDFs). Since our 
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calculation does not take into account NLO QCD effects, we consistently resort to LO PDFs, 
using the LHAPDF implementation of the central MSTW2008LO PDF set [TO2]. From there we 
infer the value 

o£°(M z ) = 0.1394 (3.8) 

for the strong coupling constant. We identify the QCD factorisation scale fip and the renormal- 
isation scale /Ur choosing 

MF = Mr = M z . (3.9) 

Note that the choice of the scales ^f,r as well as the actual value for the strong coupling a s plays 
a minor role for our numerical analysis of EW radiative corrections in Section 13.2.21 We focus on 
the relative importance of the NLO EW corrections considering the ratio cr^yp /c LO from which 
the a s - and the scale dependence drop out. 

For the jet-reconstruction we use the anti-fcx clustering algorithm [103] with separation 
parameter R = 0.4. For our scenario with exactly two partons and one potential photon 
in the final state this simply amounts to recombining the photon with a parton a if R a ^ = 
\J {Va — y^y) 2 + < R- Here y = ^ ln[(E +p\ J )/{E — pl)] is the particle's rapidity with E 
denoting its energy and j>l its three-momentum component along the beam axis, and t/> a7 is 
the azimuthal angle between the the photon and the parton a in the plane transverse to the 
beam axis. In case of recombination, the resulting photon-parton jet is subjected to the cut 
z 7 = E 7 / (E 7 + E a ) < 0.7 in order to distinguish between Z + 2 jets and Z + 1 jet + 7 production 
as explained in Section [3.1.31 After a possible recombination, we require two hard jets with 

PTjct >25GeV, M<4.5 (3.10) 

for the final event. 



3.2.2 Results 

In this section we present results for the total cross section and various differential distributions 
using the numerical input parameters and acceptance cuts introduced above. The total cross 
section and its composition at LO for the 8TeV LHC is shown in Table Q] where the absolute and 
relative contributions of the partonic channels are listed. In the lower part of Tabled] we provide 
the contribution to the total cross section of partonic processes with external gluons (gluonic) and 
of the four-quark processes (four-quark). We find the total cross section dominated by processes 
with external gluons, in particular by the quark-gluon induced processes (69%). We also provide 
the NLO cross section and the relative EW corrections for the gluonic channels in the last two 
columns of Table [TJ For our set of cuts, they range between —1.0% and —1.4% for the different 
gluonic channels. 

In the following we present results for distributions at LO and NLO for the gluonic channels 
only. Although these channels dominate the total cross section we emphasise that there are 
certain phase-space regions where the relative importance of the four-quark processes is enhanced 
and these channels and the corresponding EW corrections should not be neglected to describe 
pp — > Z + 2 jets properly. 

For each distribution we provide two plots: the upper panels show the LO and NLO prediction 
for the differential cross section while the lower panels show the NLO result normalised to the LO 
result. In Fig. [3]we present results for the differential cross section as a function of the transverse 
momentum and the rapidity for the harder jet ji and softer jet j2, respectively. Both transverse 
momentum distributions show steep slopes over six orders of magnitude in the displayed range. 
The EW corrections lower the LO prediction, and their relative size grows in absolute value with 
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Process class 


(j LO [pb] 


^ L °M L o? [%] 


*° [Pb] 


_NLO 

~ 1 [%] 


ug -)■ ugZ, dg -)■ dgZ, 
me — )■ meZ, — > <isZ 


1324.1(2) 


68.79 


1308.8(2) 


-1.16 


mm -)• ggZ, (id-)- ggZ 


128.84(2) 


6.69 


127.56(2) 


-0.99 


EE — )- UU.Z . EE — )- (i(i^ 


87.40(2) 


4.54 


86.18(2) 


-1.37 


mm' — > mm'Z, (id' — )■ dd'Z, 
mm' — y mm'Z, dd' — y dd'Z 


88.41(2) 


4.59 






uu — > u'u'Z, dd — > d'd'Z, 
mm' — > uu'Z, dd' — > dd'Z 


87.98(2) 


4.57 


— 


— 


uu —> ddZ, dd — > uuZ, 
uu' — > dd'Z, dd' — > mm'Z 


16.566(3) 


0.86 


— 


— 


ud — > u'd'Z, ud — > u'd'Z, 
ud — > udZ, ud — > udZ 


111.74(3) 


5.81 






ud — > u'd'Z, ud —7- u'd'Z, 
ud — > udZ, ud — > udZ 


79.70(2) 


4.14 






gluonic 


1540.4(2) 


80.02 


1522.5(2) 


-1.16 


four-quark 


384.41(4) 


19.98 






sum 


1924.8(2) 


100.00 







Table 1: Composition of the LO cross section for pp — > Z + 2 jets at the 8TeV LHC. In the first 
column the partonic processes are listed, where u, u' denote the up-type quarks u, c and d, d' the 
down-type quarks d, s, b. The second column provides the corresponding cross section where the 
numbers in parenthesis give the integration error on the last digit. The third column contains the 
relative contribution to the total cross section in percent. In the fourth column we list the NLO 
EW cross section for the gluonic channels and in the last column the relative EW corrections. 

increasing transverse momentum due to the well-known EW Sudakov logarithms [104H107] . For 
transverse momenta of the softer jet j>tj 2 — ITeV the impact of EW corrections amounts up 
to —10% relative to the Born approximation, while for transverse momenta of the harder jet 
PT.ji — 1 TeV the EW effects are of the order of —15%. Since, the rapidity distributions are not 
sensitive to Sudakov logarithms, the corresponding EW corrections are flat and around —1% as 
for the total cross section. 

The differential cross section as a function of the di-jet invariant mass and as a function of 
the transverse momentum of the Z boson is shown on the left-hand side in Fig. |H For both 
distributions we find the expected dependence on the Sudakov logarithms, although the sensi- 
tivity in the di-jet invariant-mass distribution is less pronounced than in the pT,z-distribution. 
For Mjj ~ 500 GeV the EW corrections are of the order of —2%; they amount up to —4% for 
Mjj~2TeV. The transverse momentum distribution of the Z boson receives large corrections, 
from —15% for pt,z — 500 GeV to —25% for pr,z — ITeV. On the upper right-hand side in 
Fig. H] we present the differential distribution of the relative azimuthal angle c/>jj between the two 
jets. The (^-distribution shows that the two jets are preferably back-to-back in the transverse 
plane and that the EW corrections lower the differential cross section by 1—1.5%. They induce 



18 





10 5 




10' 


> 

03 


in 3 

1U 


o 






10- 


-d | f 


10 1 




10° 




io- 1 




1 




-4 -3 -2 -1 



2 3 



2 3 4 -4 -3 -2 -1 

Figure 3: Distributions of the transverse momentum and the rapidity of the harder jet ji and 
the softer jet }2 at the 8TeV LHC at LO (blue, dashed) and NLO (red, solid). The lower panels 
show the ratio of the NLO distribution over the LO distribution. 



a shape change at the permille level relative to the LO approximation. The rapidity distribution 
of the Z boson is depicted in the lower right-hand side of Fig. |U In the central region \yz\ < 2, 
where most of the Z bosons are produced, the EW corrections lower the LO cross section by 
1—1.5% while for \yz\ > 2 their effect drops to the permille level. 

4 Conclusions 

The full exploitation of the Large Hadron Collider relies on precise theoretical predictions. To 
this end QCD and electroweak next-to-leading order corrections have to be calculated for many 
processes involving many particles in the final state. This requires efficient and reliable automatic 
tools. 

In this paper we have presented Recola, a Fortran90 code for the REcursive Computation 
of One-Loop Amplitudes. It uses methods based on Dyson-Schwinger equations to calculate the 
coefficients of all tensor integrals appearing in a one-loop amplitude recursively. The tensor 
integrals can then be evaluated with efficient numerically stable techniques. The algorithm has 
been implemented for the full electroweak Standard Model, including counterterms and rational 
terms, but could be generalised to more complicated theories in a straightforward way. The 
implementation supports the complex-mass scheme and is thus applicable to processes involving 
intermediate unstable particles. For the treatment of colour we have developed a new recursive 
algorithm based on colour structures that naturally appear in the colour-flow representation. 
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Figure 4: Distributions of the di-jet invariant mass, the relative azimuthal angle between the two 
jets, the transverse momentum and the rapidity of the Z boson at the 8 TeV LHC at LO (blue, 
dashed) and NLO (red, solid). The lower panels show the ratio of the NLO distribution over the 
LO distribution. 

As a first application of Recola, we have calculated the electroweak corrections to the 
dominant partonic channels in Z + 2 jets production at the LHC. The results have been verified 
with an independent calculation based on Feynman-diagrammatic methods. For a typical set of 
cuts, the electroweak corrections are negative at the level of one percent, but become sizeable 
where large energy scales are relevant. However, in general and in particular for large energies 
of the jets, a meaningful prediction requires the inclusion of next-to-leading-order corrections to 
all partonic channels. This will be pursued in a forthcoming publication. 

Acknowledgements 

We are grateful to S. Pozzorini for many useful discussions and to A. Muck for help concerning 
Pole. This work was supported in part by the Swiss National Science Foundation (SNF) under 
contract 200021-126364 and by the Deutsche Forschungsgemeinschaft (DFG) under reference 
number DE 623/2-1. 

References 

[1] M. Ciccolini, A. Denner and S. Dittmaier, Phys. Rev. D 77 (2008) 013002 |arXiv:0710.4749l 
[hep-ph]]. 



20 



[2] Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, Nucl. Phys. B 425 (1994) 217 
[hep-ph/9403226| . 

[3] Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, Nucl. Phys. B 435 (1995) 59 
[hep-ph/9409265| . 



[4] R. Britto, F. Cachazo and B. Feng, Nucl. Phys. B 725 (2005) 275 |hep-th/0412"i~03] . 



[5] G. Ossola, C. G. Papadopoulos and R. Pittau, Nucl. Phys. B 763 (2007) 147 
[hep-ph/0609007| . 

[6] G. Ossola, C. G. Papadopoulos and R. Pittau, JHEP 0707 (2007) 085 |arXiv:0704.127T1 
[hep-ph]]. 

[7] R. K. Ellis, W. T. Giele and Z. Kunszt, JHEP 0803 (2008) 003 [arXiv: 0708 .23981 [hep-ph]]. 

[8] W. T. Giele, Z. Kunszt and K. Melnikov, JHEP 0804 (2008) 049 |arXiv: 080 1.22371 [hep-ph]]. 

[9] R. K. Ellis, W. T. Giele, Z. Kunszt and K. Melnikov, Nucl. Phys. B 822 (2009) 270 
|arXiv:0806.3467l [hep-ph]]. 

[10] G. Passarino and M. J. G. Veltman, Nucl. Phys. B 160 (1979) 151. 

[11] C. F. Berger, Z. Bern, L. J. Dixon, F. Febres Cordero, D. Forde, H. Ita, D. A. Kosower and 
D. Maitre, Phys. Rev. D 78 (2008) 036003 |arXiv:0803.4180l [hep-ph]]. 

[12] W. T. Giele and G. Zanderighi, JHEP 0806 (2008) 038 [arXiv:0805.2l"52l [hep-ph]]. 

[13] A. Lazopoulos. larXiv:0812.2998l [hep-ph], 

[14] W. Giele, Z. Kunszt and J. Winter, Nucl. Phys. B 840 (2010) 214 |arXiv:0911.1962l [hep-ph]]. 

[15] S. Badger, B. Biedermann and P. Uwer, Comput. Phys. Commun. 182 (2011) 1674 
|arXiv:1011.2900l [hep-ph]]. 

[16] V. Hirschi, R. Frederix, S. Frixione, M. V. Garzelli, F. Maltoni and R. Pittau, JHEP 1105 
(2011) 044 jarXiv:1103.062T1 [hep-ph]]. 

[17] G. Bevilacqua, M. Czakon, M. V. Garzelli, A. van Hameren, A. Kardos, C. G. Papadopoulos, 
R. Pittau and M. Worek, larXiv:1110.1499l [hep-ph]. 

[18] G. Cullen, N. Greiner, G. Heinrich, G. Luisoni, P. Mastrolia, G. Ossola, T. Reiter and 
F. Tramontano, Eur. Phys. J. C 72 (2012) 1889 |arXiv:1111.2034l [hep-ph]]. 

[19] C. F. Berger, Z. Bern, L. J. Dixon, F. Febres Cordero, D. Forde, T. Gleisberg, H. Ita and 
D. A. Kosower et al, Phys. Rev. D 80 (2009) 074036 [arXiv:0907.1984l [hep-ph]]. 

[20] R. K. Ellis, K. Melnikov and G. Zanderighi, Phys. Rev. D 80 (2009) 094002 |arXiv: 0906 .14451 
[hep-ph]]. 

[21] C. F. Berger, Z. Bern, L. J. Dixon, F. Febres Cordero, D. Forde, T. Gleisberg, H. Ita and 
D. A. Kosower et al, Phys. Rev. Lett. 102 (2009) 222001 [arXiv: 0902 .27601 [hep-ph]]. 

[22] G. Bevilacqua, M. Czakon, C. G. Papadopoulos, R. Pittau and M. Worek, JHEP 0909 
(2009) 109 [arXiv:0907.4723l [hep-ph]]. 



21 



[23] K. Melnikov and G. Zanderighi, Phys. Rev. D 81 (2010) 074025 |arXiv:0910.367T1 [hep-ph]]. 

[24] G. Bevilacqua, M. Czakon, C. G. Papadopoulos and M. Worek, Phys. Rev. Lett. 104 (2010) 
162002 [arXiv: 1002.40091 [hep-ph]]. 

[25] C. F. Berger, Z. Bern, L. J. Dixon, F. Febres Cordero, D. Forde, T. Gleisberg, H. Ita and 
D. A. Kosower et al, Phys. Rev. D 82 (2010) 074002 [arXiv:1004.1659l [hep-ph]]. 

[26] T. Melia, K. Melnikov, R. Rontsch and G. Zanderighi, JHEP 1012 (2010) 053 
|arXiv:1007.53l3l [hep-ph]]. 

[27] C. F. Berger, Z. Bern, L. J. Dixon, F. Febres Cordero, D. Forde, T. Gleisberg, H. Ita and 
D. A. Kosower et al, Phys. Rev. Lett. 106 (2011) 092001 [arXiv: 1009.23381 [hep-ph]]. 

[28] G. Bevilacqua, M. Czakon, A. van Hameren, C. G. Papadopoulos and M. Worek, JHEP 
1102 (2011) 083 [arXiv: 1012.42301 [hep-ph]]. 

[29] T. Melia, K. Melnikov, R. Rontsch and G. Zanderighi, Phys. Rev. D 83 (2011) 114043 
|arXiv:1104.2327l [hep-ph]]. 

[30] R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, R. Pittau and P. Torrielli, JHEP 1109 
(2011) 061 }arXiv:1106.60l9l [hep-ph]]. 

[31] H. Ita, Z. Bern, L. J. Dixon, F. Febres Cordero, D. A. Kosower and D. Maitre, Phys. Rev. 
D 85 (2012) 031501 |arXiv:1108.2229l [hep-ph]]. 

[32] G. Bevilacqua, M. Czakon, C. G. Papadopoulos and M. Worek, Phys. Rev. D 84 (2011) 
114017 [arXiv:1108.285T1 [hep-ph]]. 

[33] Z. Bern, G. Diana, L. J. Dixon, F. Febres Cordero, S. Hoche, D. A. Kosower, H. Ita and 
D. Maitre et al, Phys. Rev. Lett. 109 (2012) 042001 [arXiv:1112.3940l [hep-ph]]. 

[34] N. Greiner, G. Heinrich, P. Mastrolia, G. Ossola, T. Reiter and F. Tramontano, Phys. Lett. 
B 713 (2012) 277 [arXiv: 1202.60041 [hep-ph]]. 

[35] G. Bevilacqua and M. Worek, JHEP 1207 (2012) 111 |arXiv:1206.3064l [hep-ph]]. 

[36] S. Badger, B. Biedermann, P. Uwer and V. Yundin, arXiv:1209.0098 [hep-ph]. 

[37] J. Alcaraz Maestre et al. [SM AND NLO MULTILEG and SM MC Working Groups Col- 
laboration] , [^Xml203l803] [hep-ph]. 



[38] D. E. Soper, Phys. Rev. D 62 (2000) 014009 hep-ph/9910292 



[39] S. Becker, C. Reuschle and S. Weinzierl, JHEP 1012 (2010) 013 [arXiv: 1010.41871 [hep-ph]]. 

[40] S. Becker, D. Goetz, C. Reuschle, C. Schwan and S. Weinzierl, Phys. Rev. Lett. 108 (2012) 
032005 [arXiv: 1111. 17331 [hep-ph] ] . 

[41] G. Ossola, C. G. Papadopoulos and R. Pittau, JHEP 0803 (2008) 042 |arXiv:0711.3596l 
[hep-ph]]. 

[42] P. Mastrolia, G. Ossola, T. Reiter and F. Tramontano, JHEP 1008 (2010) 080 
|arXiv:1006.0710l [hep-ph]]. 



22 



[43] G. Heinrich, G. Ossola, T. Reiter and F. Tramontano, JHEP 1010 (2010) 105 
|arXiv:1008.244T1 [hep-ph]]. 



[44 
[45 

[46 
[47 

[48 
[49 

[50 

[51 

[52 

[53 

[54 

[55 

[56 
[57 

[58 
[59 
[60 
[61 

[62 
[63 
[64 



R. Pittau, Comput. Phys. Commun. 181 (2010) 1941 |arXiv:1006.3773l [hep-ph]]. 

A. Ferroglia, M. Passera, G. Passarino and S. Uccirati, Nucl. Phys. B 650 (2003) 162 
|hep-ph/0209219] . 

A. Denner and S. Dittmaier, Nucl. Phys. B 658 (2003) 175 [hep-ph/0212259] . 

T. Binoth, J. P. Guillet, G. Heinrich, E. Pilon and C. Schubert, JHEP 0510 (2005) 015 
[hep-ph/0504267l. 



A. Denner and S. Dittmaier, Nucl. Phys. B 734 (2006) 62 |hep-ph/050914l|. 

A. Bredenstein, A. Denner, S. Dittmaier and M. M. Weber, JHEP 0702 (2007) 080 
[hep-ph/06lT234] . 

A. Denner, S. Dittmaier, T. Kasprzik and A. Muck, JHEP 0908 (2009) 075 |arXiv:0906.T656l 
[hep-ph]]. 

A. Denner, S. Dittmaier, T. Kasprzik and A. Muck, JHEP 1106 (2011) 069 |arXiv:1103.09l4l 
[hep-ph]]. 

A. Denner, S. Dittmaier, S. Kallweit and A. Muck, JHEP 1203 (2012) 075 |arXiv:1112.5l"42l 
[hep-ph]]. 

A. Bredenstein, A. Denner, S. Dittmaier and S. Pozzorini, Phys. Rev. Lett. 103 (2009) 
012002 |arXiv:0905.0lT0l [hep-ph]]. 

A. Denner, S. Dittmaier, S. Kallweit and S. Pozzorini, Phys. Rev. Lett. 106 (2011) 052001 
|arXiv:1012.3975l [hep-ph]]. 

F. Campanario, C. Englert, M. Rauch and D. Zeppenfeld, Phys. Lett. B 704 (2011) 515 
[arXiv: 1106.407191 [hep-ph]]. 

L. Reina and T. Schutzmeier, JHEP 1209 (2012) 119 [arXiv: 11 10.44381 [hep-ph]]. 

F. Cascioli, P. Maierhofer and S. Pozzorini, Phys. Rev. Lett. 108 (2012) 111601 
|arXiv:1111.5206l [hep-ph]] . 

A. van Hameren, JHEP 0907 (2009) 088 |arXiv:0905. 10051 [hep-ph]]. 

F. A. Berends and W. T. Giele, Nucl. Phys. B 306 (1988) 759. 

J. M. Campbell and R. K. Ellis, Phys. Rev. D 65 (2002) 113007 [hep-ph/0202T76] . 

J. M. Campbell , R. K. Ellis and D. L. Rainwater, Phys. Rev. D 68 (2003) 094021 
[hep-ph/03081951 . 

C. Oleari and D. Zeppenfeld, Phys. Rev. D 69 (2004) 093004 jhep-ph/0310T56] . 

F. J. Dyson, Phys. Rev. 75 (1949) 1736. 

J. S. Schwinger, Proc. Nat. Acad. Sci. 37 (1951) 452. 



23 



[65] J. S. Schwinger, Proc. Nat. Acad. Sci. 37 (1951) 455. 

[66] A. Kanaki and C. G. Papadopoulos, Comput. Phys. Commun. 132 (2000) 306 
[hep-ph/0002082| . 

[67] C. G. Papadopoulos and M. Worek, Eur. Phys. J. C 50 (2007) 843 [hep-ph/0512l"50] . 

[68] A. Cafarella, C. G. Papadopoulos and M. Worek, Comput. Phys. Commun. 180 (2009) 1941 
|arXiv:0710.2427l [hep-ph]]. 

[69] K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B 313 (1989) 560. 

[70] F. Caravaglios and M. Moretti, Phys. Lett. B 358 (1995) 332 [hep-ph/9507237] . 

[71] G. Ossola, C. G. Papadopoulos and R. Pittau, JHEP 0805 (2008) 004 |arXiv:0802. 18761 
[hep-ph]]. 

[72] P. Draggiotis, M. V. Garzelli, C. G. Papadopoulos and R. Pittau, JHEP 0904 (2009) 072 
|arXiv:0903.0356l [hep-ph]]. 

[73] M. V. Garzelli, I. Malamos and R. Pittau, JHEP 1001 (2010) 040 [Erratum-ibid. 1010 
(2010) 097] |arXiv:0910.3l"30l [hep-ph]]. 

[74] H. -S. Shao, Y. -J. Zhang and K. -T. Chao, JHEP 1109 (2011) 048 [ arXiv: 1106.. 50301 [hep- 
ph]]. 

[75] A. Denner, Fortsch. Phys. 41 (1993) 307 [arXiv: 0709 .10751 [hep-ph]]. 

[76] A. Denner, S. Dittmaier, M. Roth and D. Wackeroth, Nucl. Phys. B 560 (1999) 33 
[hep-ph/9904472| . 

[77] A. Denner, S. Dittmaier, M. Ro th and L. H. W ieders, Nucl. Phys. B 724 (2005) 247 
[Erratum-ibid. B 854 (2012) 504] [hep-ph/0505042] . 

[78] Z. Bern and D. A. Kosower, Nucl. Phys. B 362 (1991) 389. 

[79] G. 't Hooft, Nucl. Phys. B 72 (1974) 461. 

[80] A. Kanaki and C. G. Papadopoulos, | |hep-ph/0012004| 

[81] F. Maltoni, K. Paul, T. Stelzer and S. Willenbrock, Phys. Rev. D 67 (2003) 014026 
[hep-ph/02092711 - 

[82] A. Denner and S. Dittmaier, Nucl. Phys. Proc. Suppl. 160 (2006) 22 jhep-ph/0605312| . 

[83] D. Y. Bardin, A. Leike, T. Riemann and M. Sachwitz, Phys. Lett. B 206 (1988) 539. 

[84] W. Beenakker and A. Denner, Nucl. Phys. B 338 (1990) 349. 

[85] A. Denner, U. Nierste and R. Scharf, Nucl. Phys. B 367 (1991) 637. 



[86] A. Denner and S. Dittmaier, Nucl. Phys. B 844 (2011) 199 [arXiv: 1005.2076 [hep-ph]]. 



[87] S. Dittmaier and M. Kramer, Phys. Rev. D 65 (2002) 073007 hep-ph/0109062]. 



S. C atani and M. H. S eymour, Nucl. Phys. B 485 (1997) 291 [Erratum-ibid. B 510 (1998) 
503] [hep-ph/9605"323] . 

24 



[89] A. Denner, S. Dittmaier, T. Gehrmann and C. Kurz, Nucl. Phys. B836 (2010) 37 
|arXiv:1003.0986l [hep-ph]]. 



[90 
[91 
[92 
[93 
[94 
[95 
[96 
[97 
[98 
[99 



E. W. N. Glover and A. G. Morgan, Z. Phys. C 62 (1994) 311. 

D. Buskulic et al. [ALEPH Collaboration], Z. Phys. C 69 (1996) 365. 
A. Denner, S. Dittmaier and L. Hofer, in preparation. 

T. Motz, PhD thesis, Zurich 2011. 

T. Hahn, Comput. Phys. Commun. 140 (2001) 418 |hep-ph/0012260] . 
T. Hahn and C. Schappacher, Comput. Phys. Commun. 143 (2002) 54 [hep-ph/0105349] . 
T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118 (1999) 153 [hep-ph/9807565| . 
S. Dittmaier, Phys. Rev. D 59 (1999) 016007 [hep-ph/9805445] . 

E. Accomando, A. Denner and C. Meier, Eur. Phys. J. C 47 (2006) 125 [hep-ph/0509234] . 
S. Dittmaier and M. Roth, Nucl. Phys. B 642 (2002) 307 [hep-ph/0206070] . 



[100] J. Beringer et al. [Particle Data Group Collaboration], Phys. Rev. D 86 (2012) 010001. 

[101] [Tevatron Electroweak Working Group and CDF and DO Collaborations]. larXiv: 1 107. 5255 
[hep-ex] . 

[102] A. D. Martin et al, Eur. Phys. J. C 63, (2009) 189 [arXiv: 090 1.00021 [hep-ph]]. 

[103] M. Cacciari, G. P. Salam and G. Soyez, JHEP 0804 (2008) 063 [arXiv:0802.lT89l [hep-ph]]. 

[104] P. Ciafaloni and D. Comelli, Phys. Lett. B 446 (1999) 278 l hep-ph/9809321] . 



[105] J. H. Kiihn and A. A. Penin, hep-ph/9906545 



[106] V. S. Fadin, L. N. Lipatov, A. D. Martin and M. Melles, Phys. Rev. D 61 (2000) 094002 
[hep-ph/9910338| . 

[107] A. Denner and S. Pozzorini, Eur. Phys. J. C 18 (2001) 461 [hep-ph/001020l] . 



25 



